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Abstract 

The purpose of this paper is to present a new mathematical model 
for the dynamics of thin Cosserat elastic plates. Our approach, which 
is based on a generalization of the classical Reissner-Mindlin plate the- 
ory, takes into account the transverse variation of microrotation and cor- 
responding microintertia of the the elastic plates. The model assumes 
polynomial approximations over the plate thickness of asymmetric stress, 
couple stress, displacement, and microrotation, which are consistent with 
the elastic equilibrium, boundary conditions and the constitutive rela- 
tionships. Based on the generalized Hellinger-Prange-Reissner variational 
principle for the dynamics and strain-displacement relation we obtain the 
complete dynamic theory of Cosserat plate. 
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1 Introduction 

This paper is straightforward extension of the static theory of Cosserat plates 
[T7] for the dynamic case. In order to describe dynamics of elastic plates with 
microstructure that possess grains, particles, fibers, and cellular structures A. 
C. Eringen (1967) was the first to propose a theory of plates in the framework 
of Cosserat (micropolar) Elasticity [3] . His theory is based on a direct technique 
of integration of the Cosserat Elasticity and assumes no variation of micropolar 
rotations in the thickness direction. Eringen's plate theory in the current form 
does not produce the Reissner-Mindlin plate equations for zero microrotations. 
In this paper we propose to use the Reissner-Mindlin's plate theory as a foun- 
dation for the modeling of dynamics Cosserat elastic plates. Our approach, in 
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addition to the transverse shear deformation, takes into account the second or- 
der approximation of couple stresses and the variation of micropolar rotations 
in the thickness direction and the corresponding inertia characteristics 

2 Micropolar (Cosserat) Linear Elasticity 
2.1 Fundamental Equations 

Before proceeding some notation convention should be explained. We use the 
usual summation conventions and all expressions that contain Latin letters as 
subindices are understood to take values in the set {1, 2, 3}. When Greek letters 
appear as subindices then it will be assumed that they can take the values 1 or 
2. 

The Cosserat elasticity equilibrium equations without body forces represent 
the balance of linear and angular momentums of micropolar elasticity and have 
the following form [3J: 

diver = pii, (1) 
£ • er +div/x = Jip, (2) 

where the quantity er = {crji} is the stress tensor, n = {^ji} the couple stress 
tensor, u and cp the displacement and rotation vectors, p and J ={Ji} the 
material density and the rotatory inertia characteristics, and £ = {£ijk} is 
the Levi-Civita tensor, where e^fc equals 1 or -1 according as (i,j,k) is an 
even or odd permutation of 1,2,3 and zero otherwise, and £ • er = {sijk&jk} -The 
constitutive equation can be written in Nowacki's form |12j : 

er = (^ + a)7 + ( M -a)7 T + A(tr 7 )l, (3) 
A* = (7 + e)x+(7-e)x T + /?(trx)l, (4) 

which we consider with the strain-displacement and torsion-rotation relations 

7 = (Vu) T +£ • if and x = Vy?, (5) 

where quantities 7 and x 1 are the micropolar strain and torsion tensors, 1 the 
identity tensor, fi, A are the symmetric and 0, ,7,e, a the asymmetric Cosserat 
elasticity constants 

The constitutive equations can be written in the reversible form: 



7 = (jjf + a')tr+(jj/ -a')a- T + A'(trcr)l, (6) 
X = (7' + e')M + (Y-0^ T +/3(tr/x)l. (7) 

where M ' = a' = ± 7' = e> = A' = and = 

Neff 1 101 uses the notation fi c for elastic parameter a. 
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We consider a Cosserat elastic body Bq. In this case the equilibrium equa- 
tions (HJ) - (0) with constitutive formulas ([3]) - (J4|) and kinematics formulas ([5]) 
should be accompanied by the following mixed boundary 

u = u , <p = ip Q on Qi = 3B a \dB a , (8) 
er n = cr ■ n = <x , n n = n ■ n = /x Q on Q 2 = dB a , (9) 

and initial conditions 

u(x,o) = U , (p(x, o) = * c in B , (10) 
u(x,o) = U , ¥»(x,o) = * c on B Q , (11) 

where u Q , ip a are prescribed on Q%, <j Q and /x on Q 2 , and n denotes the outward 
unit normal vector to OBq. 

2.2 Cosserat Elastic Energy 

The strain stored energy Uc of the body Bq is defined by the integral [T2] : 

U C = [ W{-r, X }dv, (12) 
Jb„ 

where 



H + a a — a A 

W{~1,X\ = — — 7«7ij + — — lijlji + -jlkklnn (13) 

7 + e 7 — e /3 



then the constitutive relations © - Q can be written in the form: 

a = C a [W\=V y W&nd f i = C fl [W\ = V x W. (14) 
The function W is non-negative if and only if [T2] 

H > 0, 3A + 2/i > 0, 

7 > 0, 3/3 + 2 7 >0, (15) 

a > 0, [i + a > 0, 

e > 0, 7 + 0O. 

For future convenience, we present the stress energy 



Uk = / ${cr,/x}du, 
J Bo 



where 



fi' + a' [j! — a' X' 

5 > i cr I M/ = 2 a ij a ij~\ 2^ <J ii cr ji + y ff tfe <T ra 



f + e'.. , 7 / -£ / /?' 

2 /'..,/';- + "2" 



(16) 
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The reversible constitutive relation © - ([7]) can be also written in form 

der ' * = K > ] = ^ 



7 = K 7 [<r] = — , x = K x [/i] = — . (17) 



The total internal work done by the stresses er and /j. over the strains 7 and 
X for the body £? [H] is 

£7= / [<r- 7 + (18) 

J B 

and 

f/ = U K = U c 
provided the constitutive relations ((3j) - (j4|) hold. 

2.3 The Generalized Hellinger-Prange -Reissner (HPR) 
Principle 

The HPR principle [5] in the case of Cosserat elasticity states, that for any set A 
of all admissible states s = [u, tp, 7, x, tr, lA that satisfy the strain-displacement, 
torsion-rotation relations ([5]) and the initial condition, the zero variation 

SO (s) = 

of the functional 



6(b) = U K 



/ [er • 7 — pii • u + fi ■ x — J ¥> ' <p] dv (19) 
J Bo 

+ [cr n • (u - u ) +fi n (ip- tp a )] da+ [a- ■ u + m Q • ip] da 

at s GA is equivalent of s to be a solution of the system of equilibrium equations 

CJ - (EI) , constitutive relations ([6]) - (O, which satisfies the mixed boundary 
conditions (J8j> - ([9]). The proof is similar to the proof for HPR principle for 
classic linear elasticity [5]. 



3 The Cosserat Plate Assumptions 

In this section we formulate our stress, couple stress and kinematic assumptions 
of the Cosserat plate. The set of points P = {Y x [-h/2, h/2}} I) TUB forms 
the entire surface of the plate and {r„ x [— h/2, h/2]} is the lateral part of the 
boundary where displacements and microrotations are prescribed. The notation 
T a = T\T U of the remainder we use to describe the lateral part of the boundary 
edge {r^ x [—h/2, h/2]} where stress and couple stress are prescribed. We also 
use notation Pq for the middle plane internal domain of the plate. 
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In our case we consider the vertical load and pure twisting momentum 
boundary conditions at the top and bottom of the plate, which can be writ- 
ten in the form: 



Cr 3 3U=V2 = <y t (x 1 ,X 2 ,t), 33\x 3 = -h/2 = <T b (x!, X 2 , t) , (20) 

O'30\x 3 =h/2 = 0, (T3p\ X3 =-h/2 = (21) 

t*3a\x 3 =h/2 = ^( x l> x 2,t), M33|x 3 =-h/2 = f^(xi, X 2 ,t), (22) 

Pzp\x 3 =h/2 = = 0, M 3/5 U 3 = -/ i /2 = 0, (23) 

where (x\,x 2 ) £ Pq- 



3.1 Stress, Couple Stress and Kinematics Assumptions 

Our approach, which is in the spirit of the Reissner's theory of plates [13] . 
assumes that the variation of stress a k i and couple stress p, kl components across 
the thickness can be represented by means of polynomials of a^.We adapt the 
expressions for the stress and couple-stress components in the following form 



cr a /3 = n a p(xi,X2,t) + -( 3 m al 3(xi,x 2 ,t), 


(24) 


03/3 = q l3 (x 1 ,x 2 ,t) (l - C3) i 


(25) 


CT /33 = qp(xi,x 2 ,t) (l - C3) ■ 


(26) 


CT 33 = ~ Qd-Cs^l' + O'O, 


(27) 


M Q /3 = (1 -Cl)r a fi(xi,x 2 ,t). 


(28) 


M/33 = ( 3 s* (x 1 ,x 2 ,t) + m*p(x 1 ,x 2 ,t). 


(29) 


M3/3 = 0. 


(30) 


M33 = C 3 « + *> 


(31) 



where p = a t (xi,x 2 ,t) - cr b (xi,x 2 ,t), er = \ (fT*(xi, x 2 , t) + a b (xi, x 2 , t)) , 
v(xx,x 2 ) = \ (n t {x 1 ,x 2 ) - fi b (x 1) x 2 )) andt(xi,x 2 ) = \ (//(xi,x 2 ) + ^ b {x u x 2 )) 
satisfy the boundary condition requirements. We note that expression (|2"Tj) is 
identical to the expression of 033 given in [T3] in the case of a b = 0. 

We also assume displacements u a are also distributed linearly over the thick- 
ness of the plate [3] and that w 3 does not vary over the thickness of the plate, 
i.e. 



U a = U a (x 1 ,X 2 ,t) - -( 3 V a (xi, X 2 ,t), 

u 3 = w(xi,x 2 ,t), 



(32) 
(33) 
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where the terms V a (xi,X2 7 1) represent the rotations in middle plane. The vari- 
ation of microrotation with respect to x 3 be represented by means of the second 
and third order polynomials [T7] : 



<p a = Q° a (x u x2,t)(l-Cl), (34) 

^ 3 = Q%(x 1 ,X2,t) + C 3 (l-~CljQ 3 (xi t X2,t). (35) 

where C 3 = f^3, a, (3 e {1,2}. 

We also assume that initial condition can be presented in the similar form, 
so they can be reduced to the form 

U a (x u x 2 ,0) = U°(x u x 2 ),V a (x 1 ,X2,0) = V°{x 1 ,X2), (36) 

w(xi,X2,0) — w (xi,x 2 ) (37) 



and 



U a (xi,X2,0) 

w(xi,x 2 ,0) 



e° a (x 1 ,x 2 ,o) 

Q° 3 ( Xl ,X2,0) 



U°(xx,X2),V a (Xi,X2,0) = V°(xi,x 2 ), 
w°(xi,x 2 ) 



Q° oa (x 1 ,x 2 ), 

6° 3 (xi,a;2) and 9 3 (xi,x 2 ,i) = 6 o3 (a;i, x 2 ). 



(38) 



(39) 



e° a (x 1 ,x 2 ,o) 

G° 3 (x 1 ,x 2 ,0) 



Q° oa (xi,x 2 ), 



(40) 



9 o3 (a;i,a;2) and Q 3 (x 1 ,X2,t) = 6 o3 (a;i, x 2 ). 



4 Specification of HPR Variational Principle for 
the Dynamics of Cosserat Plates 

The HPR variational principle for a Cosserat plate is most appropriately ex- 
pressed in terms of corresponding integrands calculated across the whole thick- 
ness. We also introduce the weighted characteristics of displacements, micro- 
rotations, strains and stresses of the plate, which will be used to produce the 
explicit forms of these integrands. 

4.1 The Cosserat plate stress energy density 

We define the plate stress energy density by the formula; 

^ = \ dCs- (41) 



G 



Taking into account the stress and couple stress assumptions (|2~4"|) - (j3"Tj) and 
by the integrating 4> {a, jx} with respect (3 in [—1,1] we obtain the explicit 
plate stress energy density expression in the form [I7| : 



d>(S) 



A + ^ 



2V(3A + 2fi) 
A 

2ft/t(3A + 2/i) 
a + /i 



12 9 
— M 

12 



ft 



12 



8fta/i 
3(a — /t) 
lOfta/x 
3 



(1 - <M ( K + ^M aj} 



{QaQa + Q @Q ft) 



5 10 

?a<£ + -n 12 n 2 i + -Mi 2 m 2 i 

6 h £ 



3A 



5 V(3A + 2/x) 



Ta,a M PP 



5ft 7 (3/3 + 2 7 ) + ^ - + Wh G "7 ) ^ 



17ft(A + /t) 
"280/i(3A + 2/t) 

/7,(A + /t) 2 



(Qa,a) 



A 



7 + e 



~2^(3A + 2/i) hje 

P p . , M/? + 7) 
"2 7 (3/3 + 2 7 ) Kaat 



2/t(3A + 2/i) 
_3_ 
2ft 2 



(iVaa) cr 



f— .sv;.S', + ^(i-<W)# 2 



/3 7 



t 2 + M^Lv 2 



2 7 (3/3 + 2 7 ) 67(3/3 + 2 7 ) 
where the Cosserat stress set 

S = [M a p, QcQlai Raf3, Sp, N a p, M*] , 

where 



(42) 



(43) 



M af3 

Qa 
Rap 



h V f 1 r if h " 



2/i 



ft / At 2k * 



ft f 1 2ft 

2 y M Q| 3«C3 = y rQ <3' 



2 „i 



y Ma3<3 = K) 



(44) 



Here M n and M 2 2 are the bending moments, M 12 and M 2 i the twisting 
moments, Q a the shear forces, Q* a the transverse shear forces, R\\ and -R22 the 
micropolar bending moments, R12 and -R21 the micropolar twisting moments, 
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5* the micropolar couple moments, all denned per unit length, Nu and N22 are 
the bending forces, N12 and N21 the twisting forces, M* the micropolar shear 
couple-stress resultants. 

Then the stress energy of the plate P 



where Po is the internal domain of the middle plane of the plate P. 

4.2 The density of the work done over the Cosserat plate 
boundary 

In the following consideration we also assume that the proposed stress, couple 
stress, and kinematic assumptions are valid for the lateral boundary of the plate 
P as well. 

We evaluate the density of the work over the boundary r u x \—h/2,h/2] 



Taking into account the stress and couple stress assumptions (p4|) - (|3Tj) and 
kinematic assumptions (|32|) - (f35|) we are able to represent Wi by the following 
expression [17] : 



Wi = s n u =M a v a + Q*w + R a n° a + s*n 3 + N a u a + m*o°, (47) 




(45) 




(46) 



where the sets S n and U are defined as 



u 



n 



[M a ,Q*,R a ,S*,N a ,M*] , 

[* a ,w,n° a ,n 3 ,u a ,tii] 



and 



M. 



S* = 



a 



SpUp, N a = N a(3 nf3,M* = Mpnp, 



In the above np is the outward unit normal vector to r u , and 




Here ^ a are the rotations of the middle plane around W the vertical 

deflection of the middle plate, ft® the microrotations in the middle plate around 
Xk axis, U a is the in-plane displacements of the middle plane along x a axis, 3 
the rate of change of the microrotation ip 3 along X3. 

We also obtain the correspondence between the weighted displacement and 
the microrotations (j48|) and the kinematic variables by applying (??) and (|35[) 
in integration of expressions (|48p : 



* Q = V a (xx,x 2 ,t),W = w(x 1 ,x 2 ,t), (49) 

0° = kte° a (x u x 2 ,t), ^3 = y e 3 (x u x 2 ,t), 

U a = U a (xi,x 2 ,t), fl° = &l(xi,x 2 ,t), 

where coefficients k\ and k^ depend on the variation of microrotations. Under 
the conditions (|35|) we have that fef = 4 and k* — s 



5" 



The density of the work over the boundary x [— /i/2, /i/2] 
h f 1 

>V 2 = - J {cr oa u a + m oa ip a ) n a d( 3 
can be presented in the form [17j 

w 2 = s -u =n OQ * Q + n o3 w + M oa n° a + M G * 3 n 3 + s , Q c/ Q + r o3 n° 3 , 

where 

S = [n OQ , n o 3 , M oa , M o3 , 

w = [* a ,w,o°,n 3) c/ Q o°] , 

M a pnp — Il OQ , R a (jnfj — M oa , 
Q> a = n o3 , S*n a = M G * 3 . (50) 
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N a pnp = S a , (5!) 
M*n a = T o3 . (52) 



Now np is the outward unit normal vector to r^, and 

h f 1 h f 1 

n o 3 = g y (< 7 o3-vo)d( 3 , M* 3 = - J {n o3 - tn 3 )d( 3 , 



-l 



= g / CToq: ^3! t °3 = ( 77 ) / ( 3 (Vos ~ C 3 v) d (s- (53) 



We are able to evaluate the work done at the top and bottom of the Cosserat 
plate by using boundary conditions (f2"TJ)) and (|2"2l 

(c o3 W3 + m 3^ o3 )n 3 (ia = (pW + vtt 3 )da. 



TUB 

Po 



4.3 The Cosserat plate internal work density 

Here we define the density of the work done by the stress and couple stress over 
the Cosserat strain field: 



h f 1 

W3= 2J ( CT -7 + /*-xMC 3 - 



(54) 



Substituting stress and couple stress assumptions (|2"4")) - (lUT]) and integrating 
expression (|54|) we obtain the following expression [17] : 



W 3 = S ■ £ =M a pe a p + Q a u a + Ql a (j% + R a pT a p + S*T3 a + N a pv a p + M*T3 a , 

(55) 

where £ is the Cosserat plate strain set of the the weighted averages of strain 
and torsion tensors 



£ = [e a p, up, w*, t 3a , t a p, v a p, r 3i , 
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Here the components of £ are 

:! /•' 

e a /3 — 



h j CalapdCa, (56) 



w « = 4 



/ 7o3 (1 " C 2 ) dCs, (57) 



< = I I ( I -C J j^:, (5S) 



T 3a = ^3X3^3, (59) 

<s = I [_ X a0 (l-C 2 )d( 3 , (60) 



= o / 7 Q /3 d C3> (61) 



3n 



1 

2 7-i 



1 

XscAs- (62) 



The components of Cosserat plate strain (|55|) - (I5^|) can also be represented 
in terms of the components of set IA by the following formulas [17j : 



e a f3 = ^f3,a + £3a/3^3, 



T 3a = ^3,a, (63) 



r° - O 
r af3 — Si /3,a> 



V a /3 = Up <a + S 3a pQ 3 , 

We call the relation ([53")) the Cosserat plate strain-displacement relation. 

4.4 The density of the kinetic energy 

Here we define the density of the kinetic energy: 

h 



w 4 2 



J (pii • u + J if ■ if) d( 3 



> 3) 

which can be presented in the form 



w 4 = ku ■ u=i v a y a + Po ww + i oa n° a n Q Q + j*n 3 n 3 + Po u a u a + j o3 o°q°, 

where 



li 



K— [loi Poi loa, J 3 , Pa ^03] 7 

and 



KU= 



i * a , Po w, i oa n° a , j;n 3 , Po t) a , i o3 n° 3 



where 



ph 3 

I = — , p Q = ph, I oa = k 3 J a h 7 J 3 = klJ 3 h 3 , I o3 = J 3 h, 
" 3 ~ 6' 4 " 32' 

5 Cosserat Plate HPR Principle 

It is natural now to reformulate HPR variational principle for the Cosserat plate 
P. Let A denote the set of all admissible states that satisfy the Cosserat plate 
strain-displacement relation (|6"3")) and let 8 be a HPR functional on A defined 

by 

©( s ) = U k-J (S ■£-KU-U-pW~vQ Q 3 )da+ J S - (U -U a ) ds+ J S n Uds, 

(64) 

for every s = \U, £ , S] G A, then 

(59 (s) = 

is equivalent to the following plate bending (A) and twisting (B) mixed prob- 
lems. 

A. The flexural motions system of equations: 

M a p, a -Qp = I o y 0) (65) 
Ql. a +P = PoW, (66) 

R a p,a + £3f3f (Qy — Qj) = Io/3^°p, (67) 

S* a , a + e 301 M 01 = j;n 3 , (68) 
with the resultant traction boundary conditions : 



M a pnf3 = n oQ , R a pnp = M oa , (69) 
Q*n Q = n o3 , S* a n a = T o3 , (70) 

at the part and he resultant displacement boundary conditions 
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* a = t oa , w = w a , o° = n° oa , n 3 = n o3 , 

at the part T u . 

The constitutive formulas: 

flg _ 12(A + M ) M 
~ 9M aa ~ ft, 3 /i(3A + 2/i) 

ii ^A 3A , > 

" |£ "' 331 / 1 3 jU (3A + 2 M ) M ^ + 5M3A + 2/,) (Q ^ j 

<9$ 3(a + M ) 3(a-/i) 



5$ _ 3(a-/i) c> | 3(q + m) Qc 



<9Q Q 10/ia/i a lOhafi 
d$ _ 3(a-/i) 3(a + /i) 
<9Q* 10/ia/x V 10/ia M 



o d<& 6(/3 + 7 ) D 



<9i? aQ 5/i 7 (3/3 + 2 7 ) ' 

3/3 g 

1 5/17(3/3 + 2 7 ) pti 2 7 (3/3 + 2 7 ) ' 

o 9$ 3(6-7) 3(7 + 6) D 



d$_ _ 3( 7 + e) c , 
T3a ~ dS* a ~ /i 3 7 e ' 
B. The extensional motions system of equations: 



with the resultant traction boundary conditions at T a : 



Napnp = S Q , 
M*n a = M * 3 , 



and the resultant displacement boundary conditions at r u : 

U a = U oa , 5I3 = ^3. 



The constitutive formulas: 
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a a 



dN aa V(3A + 2/t) ' 

; N (a + l)(a+l) - n | o \ a 0, (82) 



r° 



" 2 V(3A + 2/i) JV («+ 1 )(«+ 1 ) 2 /i(3A + 2/0 ' 
<9$ _ a + / 

W a/3 4/ia/. 
5$ 7 + e 



9$ a + u a — u , , „ 



fM*. (84) 



9M* 4/i7e Q ' 

We also represent the above constitutive relation in the compact form: 

£ = K [S] = K ■ S, 

where we call /C the compliance Cosserat plate tensor. 
Proof of the theorem. The variation of O(s) 

Se{s) = / {(/C [S] - £) ■ SS-KU ■ SU - SS£+pSW + v5Q 3 } da 

Pa 

+ I {SS - (U -U ) + So-SU} ds + I S n -SUds. 

We apply Green's theorem and integration by parts for S and 8U [5] to the 
expression: 

J S ■ 5£da = J S 5-U ds- j '{(M Q/3 , Q - Qp-) 5Vp + Q* a JW 

Po dP P 

+ (R a p,a + £3/37 (Q* — Qj) RaP, a) 5 fig 
+ (Sa,a + e 3f3j M f>-y) ^3 + N a p, a SUp 

+ ( M l a + eapjNfh) $n° 3 }da. 

Then based on the fact that 8U and S£ satisfy the Cosserat plate strain- 
displacement relation (f63|) . we obtain 
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6@(s) = J {(AC [5]- £) ■ 5S -S5£}da 

Po 

+ f{(M afl , a -Q f3 - Io^fj) + (Q* )Ct +p- Po W^j SW 

Pa 

+ + ^PiMpi ~ ^3*^3) <^3 + (jVa/J,a - Po U^ 5Up 

+ (M* >a + e 3l3l N Pl +v- Joslig) <Jfi§}da 
+ / 8S -{U -U )ds+ I (S -S n )-SUds. 

If s is a solution of the mixed problem, then 

8Q{s) =0. 

On the other hand, some extensions of the fundamental lemma of calculus 
of variations [5] together with the fact that U and £ satisfy the Cosserat plate 
strain-displacement relation (|63[) imply that S is a solution of the A and B 
mixed problems. 

6 Field Equations Governing Flexural and Ex- 
tensional Motions in terms of Kinematics Vari- 
ables 

For the future consideration we represent the constitutive relations in the fol- 
lowing fornjl 

vh 2 

M aa = D (W a , a + v^a', a ') + _ ^ P, (85) 
M a , a = y / fj^l + + 2iV 2 (-l) a+1 (fi 3 - *a',a)) (86) 

= ^G/i(a t 2 -2«f)n°, ia +^fiO )Q ,), (87) 

i? QQ = (n°, a + (1 - *) (n°, a + n°,, a ,)) + ^H^zll ti 



2 In the following formulas a subindex a' = 1 if a = 2 and a' = 2 if a = 1 
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ft" 



G/ t (4Z 6 — lf)h 3 
a = 12? ^ * 



(1 — V z ) 1 — 1/ 

^Va'a = (1 _ jy2) (t/a'.a + Ca,a' - 2JV 2 (tfa'.a + (-l)"^)) , (91) 



» Gl?ttl?-l?)h n 
M* = _*^_L_LL n o oj (g2) 



"6 

where we use the following technical constants [4], [6]: the Young's modulus 
B = ^( 3A+2 w the Poisson's ratio ^ = tttt^ — v, the shear modulus G = „,, E 
the flexural rigidity of the plate Z) = yifi^v' 1 ) > * ne characteristic length for 
torsion l t — J 2 , the characteristic length for bending h = \\J , the coupling 



number N — ./ — j— , the polar ratio VP = ft o , k\ = I and 



Remark: The values of K\ and K2 depend on the form of approximation. For 
instance, in the Mindlin's case of dynamics [H] the value of 
slightly different. 

After substitution (JS5J> - ((Ml) into (JBSJ) - ([55 ]) and ([77 ]) - ([75] ) we obtain bend- 
ing and twisting governing systems. We write system for the flexural motions 
in the form: 



L(3 x )H-F=^,xeP , (93) 



where L (d x ) = L (j-) , 



L(€)=L(CJ = 



ill 


1*12 


Liz 


Lx4 





^16 


il2 


L22 


L23 


-^24 


—Liq 





—£13 


— L23 


L33 





L35 


-^36 


—1/14 


L24 





L44 











Ll6 


—£35 





£55 


L56 


^16 





£36 





—L 56 


L&6 


*1 


*2 W 


fi 3 


^2 


]. 





and 
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F T — [ F\ F 2 F3 F4 F§ 

In the above 



P /i 3 a*! p/i 3 a* 2 nh dw k j h 3dn 3 h T h dSl ° k T h 9n ° 2 



Lll = ill (£i,£ 2 ) = k l(l + k 2& - fc 3, L 2 2 = Lu{t2,€l), L 33 = fc 4 A, 

L44 = fc 5 A - fc 6 , L55 = £55(^1, £2) = ^1 + ^£2 - ^9; L 6G = -L 55 (£ 2 , 

L12 = fcloCl^2i L 13 = fcn^i, L14 = ki 2 £ 2 ,Li e = fci3, L23 = fell^2! 

^24 = ^512^,^35 = -fci3^ 2 , L 36 = fci 3 ^ 1 ,L 56 = fcu^^, A = ^ 

/i 2 ^(i - jv 2 ) a P _ ^(1 - jv 2 ) a P 

^1 = TnTi 7— q j ^5 ~~ 



10(l-i/) 9a;i ' 10(1 -1/) &r 2 ' 

- N 2 ) 



b 6 v ; &r 2 

Here 

h _ ^-^w-'?)a-^) , t , = 2 ^ D(1 _„), 

5/i(l-iV 2 )G/ 2 (2-*) _ 10/i(l - JV 2 )GZg _ 10/iGiV 2 

~~ 3 ' ~~ 3 ' ~~ 3" ' 

D{l + v-2N 2 ) 5Gh{2N 2 -l) , 
fcio = — g -. fc u = ^-g -, fc i2 = £W 2 (l-f), 

5GftiV 2 , 5/i(l-7V 2 )G(Z 2 (2-*)-2/ 2 ) 

«13 = 3 , «14 = ^ L - 

The correspondent boundary and initial conditions are 

T(d x )H-F*=0, xer ff , (94) 

H-H = 0, xer B , (95) 



and 



H(x,0) = H 
t H(x,O) =H 



where differential operator T (<9 X ) = T (af - ) 
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T(0 = T(O = 



and 





Tn 





T u 








T 2 i 


T22 





T 24 








T31 


T32 


T33 



















T44 




















T55 


T56 
















^66 


Fi 


F 3 * 


T 4 * 


F 5 * 


T 6 * ] 





o = [ *ol *o2 Wo ^o3 ^1 0« 2 ] 



In the above 



Tn = T 1 (^,^),T 22 =T 1 ^ 2 ^ 1 ),T 1 (^,^)=Dn 1 ^ 1 + ^^l ) n 2 ^, 



T33 — 
T55 — 



5G/i , , , ^ m GL 2 (4Z 2 - L 2 )/i 3 



5Gh 



{l 2 tni {2-^ 1 + 2lln^ 2 ),T { 



l 2 p ("l^i+na^), 
5G/i 



66 



3 



-(2? b Vei+^ 2 (2-*)n 2 5 2 ), 



L>(l + i/)(l-2iV 2 ) , L>(l + z/)iV 2 



2(1 -iV 2 ) 



1 -N 2 



D(l + v)(l-2N 2 ) £(l + t /)7V 2 
T 2 i = Dvn 2 t,i H ^ ^ ni^ 2 ,T 24 = ; 77;; — ni, 



2(1 -iV 2 ) 



1 - TV 2 



I31 = 



5Gh(l~2N 2 ) m 5G/i(l-2iV 2 ) m 5G/iA^ 2 . 

-ni,T 32 = —-7 77^"2,r 36 = — _(m -n 2 ), 



6(1 -TV 2 ) 



6(1 -iV 2 ) 
5G/i 



3(1 -iV 2 ) 



T 56 - T 65 = T 2 (e 2 ,Ci), T 2 (Ci^ 2 ) = -^-(Z t 2 ni(l - *)£ 2 + (Z t 2 - 2Z 2 )n 2 £ 1 ), 



F* = 



^ 2 



10(1 - v) 

f* = -n o3 , F 4 * - -M o3 
2GZ 2 (1 - tf ) 



"iP-n 0l ,F 2 * = - 



vh 2 



« 2 p-n o2 , 



f 5 * = 



ni t-M ol ,F* 



10(1 - 1/) 

2G/ t 2 (l - tf) 



The governing system for the extensional motions is 



n 2 i - M o2 . 



L(d x )H-F=^,x G P , 



(96) 



where 





' in 


T12 


L\3 


l (€) = £(€„) = 


T 2 i 


L 22 


L 2 3 




. T 3 i 


L32 


L33 
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and 



Here 



K4 



F T = [F 1 F 2 F 3 ], 

ii T = [Ui u 2 fi§ ] , 



Post Po~ j o3— 



d(7i - dU 2 T 
, , 



f 2 i 
T31 



«i = 



+ K 2£ 2 , ^12 = K3ClC2; r i3 = 2k4$ 2; 

T , i2,T , 22 = Tn, T23 = 2k4^ 1; 

(1-7V 2 ) 



2G an ' 2 



vki_c)oq_ p,* 
2G dx 2 ' 3 



2(1 — TV 2 ) 
\-v 

N 2 , k 5 = * 



, k 2 = 2N 2 , K3 = 1 - Ki = 

' 2 (4l 2 -l 2 )(l~N 2 ) 
2l\ 



Gh " W ' 
(1 + zv -2N 2 ) 

(1 - ") ' 



The boundary and initial conditions for the extensional system has the fol- 
lowing form: 



and 



T(0 x )H - f*= 0, x er a 
h h = 0, xer„, 



H(x,0)= H° 



9 t H(x,0)= H 
where differential operator T (<9 X ) = T (gf - ) 



t(€) = t(€ a ) = 



2ii 2i2 T 13 
T 2 i T22 T23 
T33 



(97) 
(98) 



_ \ T 



(f) 
(S.) 



Fl i*2 -F3 



[ u i u o2 n ( > 3 ] . 
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In the above 

- Ehn lt , Ghn 2 ~ Ehisni Ghn 2 {l - 2N 2 ) 

1 i - ^ 1 + I^v^ 2 ' Tl2 = T^^ + — y^Yp — ^ 

~ 2N 2 Ghn 2 - _Ehvn 2 Ghm(l - 2N 2 ) 

13 ~ i-jv 2 ' 21 ~ I - + I - iv 2 ^ 2 ' 

- Ehn 2 c Gh ni ~ 2N 2 Ghn x 

2 ~ 1 - v 2 1 + 1 - N 2 ^ 1 3 ~~ 1 - AT 2 ' 

- Gl 2 (Al? - l 2 )h „ 

1 — v 1 — V 



7 Conclusion 

We proposed a new mathematical model for dynamics of Cosserat elastic plates 
based on Reissner-Mindlin's plate theory. The polynomial approximations of 
the variation of couple stress and micropolar rotations in the thickness direc- 
tion allowed us to project Cosserat 3D Elasticity dynamics equations into the 
dynamics equations in the middle plane of the plate. We generalized fir the 
dynamics Hcllinger-Prange -Reissner (HPR) principle to derive the dynamics 
equations in the middle plane and constitutive relationships for the plate. In 
terms of the kinematic variables, the total system of dynamic equations de- 
scribes the flexural (subsystem of 6 equations) and the extensional (subsystem 
of 3 equations) motions of the plate. 
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